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Abstract 

It has been recently suggested that, as a gravitational Aharonov-Bohm effect due to a grav- 
itomagnetic potential, possible effects of Chern-Simons gravity on a quantum interferometer are 
dependent on the latitude and direction of the interferometer on Earth in orbital motion around 
Sun. Continuing work initiated in the earlier publication [Okawara, Yamada and Asada, Phys. 
Rev. Lett. 109, 231101 (2012)], we perform numerical calculations of time variation in the induced 
phase shifts for nonequatorial cases. We show that the maximum phase shift at any latitude might 
occur at 6, (and 12), and 18 hours (in local time) of each day, when the normal vector to the 
interferometer is vertical, eastbound and northbound, respectively. If two identical interferometers 
were located at different latitudes, the difference between two phase shifts that are measured at 
the same local time would be 0(s'm5(p) for a small latitude difference 5(p. It might thus become 
maximally ~ 20 percents for 5(p ~ 10 degrees, for instance. 

PACS numbers: 04.25.Nx, 04.50.-h, 04.80. Cc 
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I. INTRODUCTION 



As a fundamental problem, the interplay between the quantum theory and the gravita- 
tional physics has been studied mostly by theoretical experiments Corella, Overhauser, 
and Werner (COW) 2J succeeded in a first experiment involving both the Plank constant 
h and the gravitational constant G by using a neutron interferometer. In the COW experi- 
ments, a neutron interferometer is tilted, such that a neutron beam path I is higher above 
the surface of Earth than the other path segment II, causing a gravitationally induced phase 
shift of the neutron de Broglie waves on path II relative to path I. The gravitationally in- 
duced phase shift was experimentally observed 

as 

. This experimental result means that 
the classical gravitational field (at the Newtonian order) affects a quantum particle as well 
as a classical one. In recent decades, there has been technological progress in quantum ex- 
periments including neutron interferometers and quantum optics. As a next step, current 
attempts to probe general relativistic effects in quantum mechanics focus on precision mea- 
surements of phase shifts in quantum interferometers (e.g. |4|). Hogan has recently proposed 
an ambitious idea to use quantum interferometers as an experimental probe of a quantum 
spacetime at the Planck scale [5|. Quantum experiments may play a role in probing an 
intermediate regime between general relativistic gravity and Planck scale physics. 

In addition to the above fundamental issue, current astronomical observations, such as the 
apparent accelerated expansion of the Universe, suggest a possible infrared modification to 
general relativity. For instance, dark energy is introduced to explain the observed accelerated 
expansion by means of an additional energy-momentum component in the right-hand side 
of the Einstein equation. As an alternative approach for interpreting the present accelerated 
expansion, the left-hand side of the Einstein equation, equivalently the Einstein-Hilbert 
action, could be modified in various ways (nonlinear curvature terms, higher dimensions, 
and so on). The Chern-Simons (CS) correction is one of modified gravity models. The 
CS modification is not an ad hoc extension, but it is actually motivated by both string 
theory, as a necessary anomaly-canceling term to conserve unitarity [6(, and loop quantum 
gravity (LQG), as a counter term for the anomaly 7| and recently as the emergence of the 
CS gravity when the Barbero-Immirzi parameter of LQG is promoted to a scalar field and 
the Hoist action is coupled to fermions [8]. Alexander and Yunes have recently pointed out 
that CS gravity possesses the same parametrized post-Newton (PPN) parameters as general 
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relativity, except for the inclusion of a new term, proportional to the CS coupling and the 
curl of the PPN vector potential BQ. The, have als o shown that this new correction 
might be used in classical experiments, such as the Gravity Probe B (GPB), to bound CS 
gravity (see [llj for an extensive review of CS modified gravity). Their proposal has been 
implemented by Ali-Haimoud and Chen 12[ to constrain CS gravity. Dyda, Flanagan and 
Kamionkowski [13] have constrained nondynamical CS gravity (or equivalently CS gravity 
with a canonical background cosmological scalar) by studying vacuum instabilities. 

It is interesting to study, as an attempt to probe quantum gravity, possible effects of the 
CS modified gravity on quantum experiments. Along this course, Nandi and his collaborators 
have discussed the quantum phase shift in a CS modified gravity model, where an 
isolated gravitating body was considered. Their conclusion is that the induced shifts by the 
spin of the body are too tiny to be observed. However, Earth's orbital angular momentum 
(~ 3 x 10 40 kg-m 2 s _1 ) is much larger than its spin angular momentum (~ 7x 10 33 kg-m 2 s _1 ). 
3oth of the axial vectors may play a role in CS gravity. Therefore, the present authors 
15l | have considered gravitationally interacting bodies in order to investigate the quantum 
mechanical effects of Earth's orbital angular momentum in CS gravity. It has been suggested 
that a CS modified gravity theory may predict daily and seasonal phase shifts in quantum 
interferometers, which are in principle distinct from the general relativistic effects. This 
feature can be currently used as a quantum tool to probe CS gravity. 

The numerical calculation in the previous paper [15] has been done only for the equa- 
torial case for its simplicity. In the present paper, therefore, we shall perform numerical 
calculations of possible daily and seasonal variations in quantum interference at nonequa- 
torial locations. This might be important for high-precision quantum experiments, because 
of the geographical reason. Namely, almost all quantum interferometers are located not at 
the equator of Earth but at the middle latitude region between 23°26 22" and 66°33 39", 
including most areas of the USA, many European countries and Japan. We shall argue that 
the CS latitude effects are important in experiments. If some variation were marginally 
detected in the future (presumably at a low signal-to- noise ratio), it would be difficult to 
disentangle the CS signal from other effects without taking account of the latitude effect. 
Comparing phase measurements at two (or more) latitudes would be helpful for improving 
the CS bound or distinguishing the CS signal from others. Namely, a signal-to-noise ratio 
could be increased by a combined analysis of phase measurements at different latitudes. 
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II. QUANTUM INTERFERENCE INDUCED BY CS GRAVITY 



In this section, we summarize the basics of computing the quantum phase shift induced 



by CS gravity, following Ref. 



15|- 



A. CS modified gravity 

CS gravity modifies general relativity via the addition of a correction to the Einstein- 



Hilbert action, namely [16| 







Sna EE — *— / d 4 X-fR*R, (1) 

where / is a prescribed external field (with units of area in geometrized units), R is the 
Ricci scalar, and the star stands for the dual operation, such that 

R*R = -Rap~ f &£ al3liv R 1& nv, (2) 

with e^s-y the totally-antisymmetric Levi-Civita tensor and R lMV s 1 the Riemann tensor. We 
concentrate on a nondynamical (kinematical) model of CS modified gravity, where we assume 
that / does not have dynamics (no kinetic term). We focus on a slowly changing / where 
/ is considered but no second-time derivatives of / appear. The nondynamical CS theory 
is tractable and could become a good approximation in weak fields, though there remains 
a possible evolution problem of the external field / (presumably near the central region) 
consistent with the Pontryagin constraint and recent studies su gges t that the nondynamical 



CS gravity model is an overconstrained system of equations [18L Il9j. A full dynamical st udy 
of seeking approximate solutions for rotating extended bodies has yet to be carried out [11] . 
It has been accomplished in the slow-rotation approximation |2o| . 
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The weak-field solution to the CS modified field equations in PPN gauge is given by 



Hi 



#oo = -1 + 2U- 2U 2 + 4$i + 4$ 2 + 2$ 3 + 6$ 4 + 0(6), (3) 
9oi = —Vi - \Wi + 2/(V x V) t + 0(5), (4) 
9ij = (1 + 2^7)^ + 0(4), (5) 



where U, $i, $ 2 , $3, $4, Vi, W{ are PPN potentials (e.g. [2lJ), 0(A) stands for PN remainders 
of order 0(1/ c A ) for the light speed c and the dot denotes the derivative with respect to 
x° ee ct. 
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Henceforth, we investigate the quantum mechanical effects at the linear order of the / 



term in Eq. 



. See [15] for discussions on the second (or higher) order contributions. 
Following [9|, let us consider a system of nearly spherical bodies in the standard PPN 
point-particle approximation, since we concentrate on weakly-gravitating bodies. For the 
above vector potential Vi, the CS correction to the metric becomes in the barycenter frame 
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(6) 



with m A the mass of the Ath body, r A the field point distance to the Ath body, n\ = x A /r A 
a unit vector pointing to the Ath body, v A the velocity of the Ath body, J A the spin- 
angular momentum of the Ath body, and the • and x operators are the flat-space inner and 
outer products. Note that the CS correction couples with the spin and the orbital angular 
momenta. 



B. Quantum phase shifts in a curved spacetime 

We consider a quantum interferometer that consists of a closed path C (its area S) on 
Earth, as shown by Fig [TJ 

The Hamiltonian for a quantum particle in a curved spacetime involves g^. The linear- 
order correction to the Hamiltonian by go? becomes 5H = mcgoiV 1 for a slowly moving 
particle j^. A phase difference induced by is thus expressed as ^ 

A = i / SHdt 

TYIC f 

= -jrfj-df, (7) 

where g denotes (g i, <?02> flto), m denotes the quantum particle mass, h = h/2n denotes 
Dirac's constant. By using Stokes theorem, A is rewritten in the surface integral form over 
S as 

A = ^ f(Vxg)-dS. (8) 



h 



s 



This form has the exact analogy in the Aharonov-Bohm (AB) effect, because g is the grav- 
itomagnetic potential in the gravito-electromagnetic description. The original AB effect in 



the phase shift, which was confirmed experimentally 25|. is oc § C A- dr= J S (V x A) ■ dS for 
a vector potential A in the electromagnetism. Note that the phase difference A in Eq. (JSJ) 
is caused by frame- dragging effects in rotational spacetimes and hence it does not depend 
on de Broglie wavelength A, in contrast to COW experiments. 



C. Phase shifts induced by CS gravity 

Let us substitute the CS term of Eq. (JBJ) into Eq. (JBJ) to obtain A for CS gravity. By 
using an identity e 1 -* (l/r)jki = with the Levi-Civita symbol e ijfc , one can see that the J- 
dependent part of the metric in Eq. (JBJ) always vanishes in Eq. (JSJ), whereas the v-dependent 
part makes contributions. 

Since A involves the curl operation on the surface of Earth and Earth's radius r E is 
much shorter than 1AU, the terms associated with the solar mass M & in Eq. (JSJ) are 
O(M M^ 1 r| ; lAU _3 ) ~ 10~ 9 smaller than those associated with Earth's mass M E , so that 
the terms with the solar mass (and other planetary ones) can be safely neglected. Henceforth, 
we focus on Earth's mass (also its spin and orbital angular momentum) in CS gravity. Hence, 
Eq. (JHJ) becomes [l5| 

A CS = / f^[S(v E -n E )n E -v E }.NjdS + 0(f 2 ) 
J s 

where we used r E ^> v^S* (Earth's radius is much larger than the interferometer arm length) 
and hence r = r E in the integrand. Here, in an inertial frame, v E denotes Earth's orbital 
velocity, n E stands for the unit vertical vector on the ground (at a certain latitude), Nj 
means the unit normal to the interferometer plane (See also Fig. [T|). The unit normal 
vectors n E and Nj in an inertial frame change with time as Earth rotates. The change rate 
depends on the latitude. Moreover, iVj depends also on the interferometer's direction such as 
horizontal and vertical. In contrast to COW experiments, the interferometer direction such 
as North and East does matter in CS gravity. Therefore, the factor [3(v E ■ n E )n E — v E ] ■ Nj 
in Eq. ([9]), depending on the latitude and direction, changes with Earth's spin and orbital 
motion. The directional dependence of the general relativistic phase shifts (e.g. 
in principle distinct from that of CS gravity [l^. 
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27|) is 
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The order-of-magnitude estimation of Aqq has been made in Ref. 15J. The magnitude 
of Eq. Qj is factored as 

i^-C£)(^?)(S). 

where [3(ve ■ n^n^ — ve] ■ Ni ~ 2ve- It is worthwhile to mention that the first fraction 
in the right-hand side of Eq. ( flOl) is due to the quantum mechanical physics and it is large 
enough ~ 10 24 s _1 to compensate for the factor in the second parentheses due to the CS 



ast factor in Eq. (flQj) 
2811 ) to Earth's radius. 



gravitational effect ~ fc^ 1 x 10~ 14 , where m is neutron mass. The 
is the squared ratio of the interferometer arm length (often ~ 60 cm 
In total, the magnitude of Aqq is 

'-csl -urv x (J^) (Z) (^) . (ii, 

The current bound on t he f parameter by neutron interferometry and a possible improve- 



ment have been discussed 15j . Current measurements of phase shifts in neutron interferom- 
etry do not report any anomalous (daily nor seasonal) variations with phase measurement 
accuracy at O(10~ 3 ). Current neutron interferometry, therefore, places a bound on CS 
gravity as fc^ 1 < 10°s (/ < 10 5 km), which is worse by three digits than the constraint 
fc~ 1 < lCT 3 s by the classical experiment GPB (Gravity Probe B) [9j, |29j and also LAGEOS 



30]. 



Future progress in quantum technology may improve the bound. A bound comparable 
to the GPB limit would be placed, if neutron interferometry were sufficiently improved for 
A x S^ 1 (nearly by three digits). For instance, ~ 5 meters arm length and ~ 10~ 4 phase 
measurement accuracy are preferred. 

Werner and his collaborators have already obtained the result for the measured phase 
shift with a one-sigma statistical error bar as ±0.34 mrad ~ O(10~ 4 ), where approximately, 
a total of 500 000 000 neutrons were counted in the interferograms over a period of 2 years 
(See [31] for a review of observations of Aharonov-Bohm effects by neutron interferometry). 
Furthermore, Seki and his collaborators have recently developed a multilayer cold-neutron 
interferometer and experimentally the phase measurement accuracy of 0.01 rad, where they 



used only 1.5 x 10 5 neutrons in a short time ~ 49 hours 32|. Motivated by quantum compu 



tations, for instance, Pushin and his collaborators have demonstrated experimentally how 



quantum-error-correcting codes may be used to improve experimental designs of quantum 



33|]. 



devices to achieve noise suppression in neutron interferometry 

Experimental setups usually suffer from many other seasonal variations. Lacking a signal, 
a constraint may be placed on /. In the presence of a signal, on the other hand, one would 
have to eliminate all other possible sources of seasonal variability. 



III. NUMERICAL CALCULATIONS 

A. Calculations in the laboratory frame 

The time dependence of Aqq comes from the factor [3(v E • n E ^)n E — ve] • Ni in Eq. 
09]). For computing the time variation, we take account of Earth's parameters such as the 
inclination angle of Earth's axis I E , the mean orbital angular velocity Q E , the spin rate u E , 
whereas the eccentricity of Earth's orbit has a tiny input to be ignored in this paper. 

By straightforward calculations, the above factor is rearranged as 

[3(v E ■ n E )n E - v E ] ■ Nj = {TZ^v E ) T ^3(n E0 ■ N I0 )n E0 - iV /0 ] . (12) 

Here, the superscript T denotes the transposition of the vector (and matrix), the subscript 
denotes the initial time that is chosen as the midnight on the winter solstice day and 1Z~ X 
denotes symbolically the inverse matrix of the product of rotational matrices, each of which 
corresponds to the spin of Earth, the inclination of Earth's rotation axis, its orbital rotation 
around the barycenter and the latitude angle (p, respectively. To obtain Eq. (fl2l) . we have 
used 1Z T = TZ^ 1 and a T lZb = (TZ~ l a) T b for vectors a and b. 

Equation ffl2|) seems convenient, since iV/o is fixed in the laboratory frame and the com- 
bination of [3(n E Q ■ Nj )n E o — Nj ] is thus a constant vector (as usually taken for granted 
in local experiments). On the other hand, [3(n E ■ Nj)n E — Ni] changes with time in inertial 
frames such as the barycenter frame. For discussing the time evolution of Eq. ({12"]) . it is 
sufficient to consider only the rotational matrix part. 

For treating the vector components in numerical calculations, we adopt the Cartesian 
coordinates (x, y, z) in the laboratory frame associated with the quantum interferometer, 
such that the x, y and z axes can be along the East, North and vertical upward directions, 
respectively. Hence, n E = (0,0,1). In the laboratory frame, the Cartesian components of 



S 



the term 1Z 1; ve become 



(TV Ve) x =ve ([sin 2 (f2£;t) cos(Jb) + cos 2 (f2£t)] cos{uEt) 

— [sm{yL E t) cos(VL E t)(l — cos(Je))] sin(a; E t)) , 

(TZ~ 1 ve) v =ve (— sm(Lp)[{sm 2 (Q E t) cos(Ie) + cos 2 (f2£/!;)} sin(u; E £) 
+ {sm^gi) cos(nEt)(l — cos(/b))} cos(u;e£)] 
+ cos((/?) sin(f2 E t) sin^)) , 

(JZ~ 1 ve)z =ve (cos(</9)[{sin 2 (fi£;t) cos(Je) + cos 2 (f2£t)} sm(ujEt) 
+ {sin(flEt) cos(flEt)(l — cos(Ie))} cos(cuEt)] 
+ sin(ip) sm(Q E t) sm(I E )) ■ (13) 

Let us investigate the time at which the phase shift becomes maximum each day. The 
expression of the above factor is linear in both sin(u>Et) and cos(o;£;t), because it is linear in 
the velocity ve- Therefore, it is shown that <9(AQg)/<9t = is approximately equal to 

tan(u E t) = F(Q E t, tp, I E , N I0 ). (14) 

— * 

Here, F is a lengthy function of VL E t, <p, Ie and N I0 without including uj e an d ve, and we 
have used HeCwb that implies \d (sin flEt)/ 'dt\ « \d (sin uiEt)/dt\. 

Let us consider three cases that Nj is vertical upward, eastbound and northbound. The 
extremum condition Eq. (TT4")) becomes simpler and it implies that, independent of (p, the 
maximum phase difference occurs at 6, (and 12), and 18 hours of each day, when the 
normal vector to the interferometer is vertical, eastbound and northbound, respectively. 
This result agrees with numerical ones. 



B. Equatorial cases 



For its simplicity, Ref. [15| made numerical calculations only for the equatorial case. 
Figures [2] and [3] show numerical results of time variations in the phase difference. For the 
northbound case, one-day variation is extremely small. Hence, the dashed curve for Aqq in 
Figure [2] looks like a horizontal straight line. 

Aqq on the winter solstice day is the same as that on the summer solstice day. This is 
because the angles between Earth's axis and its orbital velocity on the two days agree with 



9 



each other. On the other hand, Aqq on the spring equinox day and that on the autumn 
equinox day are the same with the opposite sign, since the relative angles on the days are 
too. 

C. Nonequatorial cases 

For the case of if = 45°, Figure H] shows variations of the phase shift on the winter 
solstice day, the vernal equinox day, the summer solstice day and the autumn equinox day. 
For two cases as Nj being along the vertical direction and along the East one, time variation 
behaviors are slightly different from the equatorial cases in Figure [21 whereas, for Nj pointing 
the North direction, they are significantly different from the equatorial case. Namely, the 
phase shift for the northbound Nj case depends more strongly on the latitude. Eqs. f|T2|) and 
(|T3l) suggest that Aqq is linear in both sin(^) and cos(^). For a small latitude difference 
Sip, therefore, the change in Aqq would be 0(Sip) ~ 0(sm5ip). 

The seasonal variation for if = 45° is plotted in Figure Like the seasonal variation 
of the weather such as the maximum and minimum temperatures, the seasonal variation 
of Aqq is more strongly dependent on the latitude than the daily variation. On all the 
four days, the change of Aqq for the northbound case is smallest. More generally, Figure 
suggests that, the one-day variation for the northbound case is always smaller than those 
for the other two cases. Therefore, the northbound case might be preferred in testing a CS 
model by quantum interferometry. 

What is a role of the CS latitude effect in experiments? Without any anomaly in phase 
shifts measured by experiments, detailed theoretical templates of the time variations would 
not be needed. The latitude effect seems less important for this case. On the other hand, 
if some variation were marginally detected in the future (conceivably at a low signal-to- 
noise ratio), it would be difficult to distinguish the CS signal from other effects without 
taking account of the latitude effect. Comparing phase measurements at two (or more) 
latitudes would be helpful for improving the CS bound or distinguishing the CS signal from 
others. Namely, a signal-to-noise ratio could be increased by a combined analysis of phase 
measurements at different latitudes. 
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IV. CONCLUSION 



For investigating possible latitude effects of Chern-Simons gravity on quantum interfer- 
ence, we performed numerical calculations of time variation in the induced phase shifts for 
nonequatorial cases. At any latitude, the maximum phase shift might occur at 6, (and 
12), and 18 hours (in local time) of each day, when the normal vector to the interferometer 
is vertical, eastbound and northbound, respectively. 

If two identical interferometers were located at different latitudes, the difference between 
phase shifts that are measured at the same local time would be 0(sin Sip) for a small latitude 
difference 5f. It might become maximally ~ 20 percents for 5<p ~ 10 degrees, for instance. 

In the northbound case, the daily variation is very small. Therefore, the northbound case 
might be used for effectively improving the statistics by combining many data points for 
each day. 

We would like to thank N. Yunes and S. Takeuchi for useful discussions. This work was 
supported in part (K.Y.) by Japan Society for the Promotion of Science, Grant-in-Aid for 
JSPS Fellows, No. 24108. 
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at some latitude 



FIG. 1: Quantum interferometer on Earth orbiting Sun. The orbital plane is chosen as the X-Y 
plane. Earth's axis and orbital velocity are denoted by Ie and ve, respectively. Top: Earth orbiting 
the barycenter. Bottom: Interferometer at a certain time and place on Earth. The latitude and the 
longitude are specified by the unit normal He, which rotates in an inertial frame around Earth's 
axis and hence its direction changes owing to the orbital motion of Earth. In an inertial frame, the 
interferometer's direction Nj also changes as Earth rotates. 
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FIG. 2: One-day variation in phase differences by CS effects on the quantum interferometer located 
on the equator of Earth (in local time). Top left: Winter solstice day. Top right: Vernal equinox 
day. Bottom left: Summer solstice day. Bottom right: Autumnal equinox day. The vertical axis 
denotes [3(ve ■ UE)nE — ve] ' Nj in Eq. Q in the units of ve = 1. We consider three cases of the 
interferometer direction. The solid, dashed and dotted curves correspond to Nj that is vertical 
upward, eastbound and northbound, respectively. 
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FIG. 3: Seasonal variation in phase differences by CS effects on the same quantum interferometer 
in Fig. [2j The winter solstice day is chosen as month and the summer solstice corresponds to 
six months. Upper: Vertical Nj (corresponding to the solid curve in Fig. [2]). Middle: Eastbound 
Nj (corresponding to the dashed curve in Fig. [2J. Bottom: Northbound Nj (corresponding to the 
dotted curve in Fig. [2]). 
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FIG. 4: One-day variation at the latitude ip = 45°, corresponding to Figure[2]for ip = 0°. Top left: 
Winter solstice day. Top right: Vernal equinox day. Bottom left: Summer solstice day. Bottom 
right: Autumnal equinox day. The solid, dashed and dotted curves correspond to Nj that is vertical 
upward, eastbound and northbound, respectively. 
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FIG. 5: Seasonal variation at the latitude tp = 45° corresponding to Figure[3j Upper: Vertical Nj 
(corresponding to the solid curve in Fig. 0]). Middle: Eastbound Nj (corresponding to the dashed 
curve in Fig. HJ). Bottom: Northbound Nj (corresponding to the dotted curve in Fig. 
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